Quantum Entanglement in the S=l/2 Spin Ladder with Ring Exchange 
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In this paper we study the concurrence and the block-block entanglement in the S — 1/2 spin lad- 
der with four-spin ring exchange by the exact diagonalization method of finite cluster of spins. The 
relationship between the global phase diagram and the ground-state entanglement is investigated. 
It is shown that the block-block entanglement of different block size and geometry manifests richer 
information of the system. We find that the extremal point of the two-site block-block entanglement 
on the rung locates a transition point exactly due to SU(4) symmetry at this point. The scaling 
behavior of the block-block entanglement is discussed. Our results suggest that the block-block 
entanglement can be used as a convenient marker of quantum phase transition in some complex 
spin systems. 
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I. INTRODUCTION 



Entanglement, as one of the most intriguing feature of 
quantum mechanics Q , has become a subject of intense 
interest in recent years. Besides being recognized as a 
kind of crucial resource of quantum computing and quan- 
tum information process [3, Q > it has also provided new 
perspectives in problems of various many-body systems. 
Particularly, the entanglement can well characterize the 
features of quantum phase transition (QPT) |_| . Many 

works ii&aiiiEiiiiiiiiiiiiiaiaiiiiiiiii 

have been devoted to understanding the relationship be- 
tween QPT and the entanglement in different systems. 
It has been observed that quantum phase transitions are 
signaled by critical behaviors of concurrence [2fJ| , a mea- 
sure of entanglement for two-qubit system, in a number of 
spin models H, IS IS El • For example, it was reported 
that the first derivative of the concurrence diverges at 
the transition point in the one-dimensional transverse 
field Ising model while the concurrence shows cusp- 
like behavior around the critical point in some 2D and 
3D spin models 0. Besides the concurrence, the block- 
block entanglement |'2l| which involves more system de- 
gree of freedom was introduced 0, 0, 0] . Especially in 
fermionic systems in which the concurrence is usually not 
applicable, the block-block entanglement can also mani- 
fest interesting properties, such as logarithmic di verg ence 
in the critical region, in a certain class of models |l3lll^| . 

However, most of the previous works were restricted 
within the models with two-body interaction, the entan- 
glement in the models with three-body or four-body in- 
teractions is less investigated and understood. In fact 
a system with multi-body interaction is important both 
in quantum information theory and condensed matter 
physics. It was pointed out that a small cluster of spins 
with three-body or four-body interactions such as the 
four-spin ring exchange could be used for quantum com- 
puting |23ll24| . Moreover, four-spin ring exchange exists 
in many physical systems and plays an important role 
in understanding the magnetism in several 2D quantum 
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FIG. 1: A sketch of spin ladder with ring exchange. 



solids such as solid 3 He [25| and Wigner Crystals |2q |. 
Therefore, it is of importance to study the properties of 
the entanglement in those spin systems with multi-body 
interactions. 

In this paper, we consider a two-legged S = 1/2 lad- 
der with additional four-spin ring exchan ge. The system 
has a very rich phase diagram |23, |28|, |2ti H3 with many 
exotic phases. We investigate the concurrence and the 
block-block entanglement in this system, and try to re- 
late them with the global phase diagram. The rest of the 
paper is organized in the following way. In section^ wc 
introduce the model Hamiltonian and its phase diagram. 
In section llTTI we study the ground-state concurrence and 
discuss its relationship with the phase diagram. In sec- 
tion IIVI we show that the two-site block-block entan- 
glement is exactly either maximal or minimal at a QPT 
point. In section [3 we show that the scaling behavior 
and some extremal point in the block-block entanglement 
can be used as marker of QPTs. In the final section, wc 
summarize our results and draw conclusions. 



II. MODEL HAMILTONIAN AND PHASE 
DIAGRAM 

The two-legged S =1/2 spin ladder with ring exchange 
(as shown in Fig. ^) is described by the following Hamil- 



tonian 



i i 



where i = 1, . . . , N/2, N is the total number of spins, 
Si.i (S2,i) is 1/2 spin operator on the upper (lower) leg 
at the ith position, and J/ (J r ) is the bilinear exchange 
constants along the legs (on the rung) and K is the cou- 
pling constant of four-spin cyclic exchange interaction P. 
Pi.i+i(Pi~i\i) rotates the four spin in the ith plaquette 
clockwise (counterclockwise), i.e. 



P 



a b 
d c 



d a 
c b 



and P 



a b 
d c 



b c 
a d 



and they can be decomposed in terms of spin operator 
involving bilinear and bi-quadratic terms, 



P + P- 



- + S a Sb + SbS c + S C S 

(s a s b ) (s c s d ) + (s a s d ) (s b s c ) 



(^SbSdJ 



(2) 



Jr 



Following the convention in Ref. [28|, we set J; 
cos 9 and K = sin 9 in the following calculation. 

Previous studies [28| suggested a rich phase diagram 
in the parameter space of 9 shown in the Fig. [3 There 
are typically six phases (regions): the rung singlet phase, 
the staggered dimmer phase, scalar chirality phase, dom- 
inant vector chirality region, dominant collinear spin re- 
gion and the ferromagnetic phase. Squares in Fig. de- 
note first-order phase transitions, circles denote second- 
order phase transitions, and the dashed line indicates a 
crossover boundary without a phase transition. 

Using the exact diagonalization method, we obtain the 
ground-state concurrence and the block-block entangle- 
ment in a spin ladder up to N = 12x2 sites with period- 
ical boundary conditions. Although the staggered dim- 
mer and scalar chirality phases are Zi symmetry break- 
ing phase with double degeneracy in the thermodynamic 
limit, the ground state is unique for most values of 9 ex- 
cept the ferromagnetic phase in a finite-size system. We 
select S z — state out of the N + 1 fold degenerate 
S = N/2 ferromagnetic states in the following calcula- 
tion. 



III. GROUND-STATE CONCURRENCE 

The entanglement between the spins at site i and site 
j can be measured by the concurrence [2(j. Let pij be 
the reduced density matrix which is obtained by trac- 
ing out all degrees of freedom of spins except that at 




FIG. 2: (color online) The ground-state concurrence of two 
spins on a rung, leg, and diagonal bond as a function of 9 in 
N — 12 x 2 spin ladder with ring exchange. The dashed line 
is the boundary of different phases. The squares on these line 
denote a first-order phase transition, the black circle denotes 
a second order transistion, while the shaded circle indicates a 
transition between two ajacent regions. 



sites i and j, and pij be the spin-reversed reduced den- 
sity matrix of p,^, i.e., p = (a y ® a y ) p* [a y eg) a y ), where 
<j y is the Pauli matrix. The concurrence C is given 
by C = max(Ai — A2 — A3 — A4,0), where {A^} are the 
square roots of the eigenvalues of the matrix pp in de- 
scending order. C — means no entanglement, while 
C = 1 the maximum entanglement such as that roots in 
Bell states. 

In Fig. |21 we show the ground-state concurrence as a 
function of 9 for a N — 12 x 2 system. In the ferro- 
magnetic phase, we can see that the concurrence on any 
two sites is the same and equals to 1/(N — 1). It van- 
ishes in the thermodynamic limit (N — > 00). In the rung 
singlet phase, we can observe a rather large concurrence 
(9 ~ 0.77r) between the two spins on the same rung. This 
fact is consistent with the picture that the ground state is 
approximated by the product state of spin singlet on the 
rungs. Similarly, the concurrence of two adjacent spins 
on the same leg is consistent with the physical picture 
of staggered singlets on the leg in the staggered dimmer 
phase. We notice that the peak of concurrence on the 
leg (~ 0.3) is much smaller than that in the rung sin- 
glet phase (~ 0.7). This is because the ground state in 
the staggered dimmer phase is two-fold degenerate in the 
thermodynamic limit. Then in a finite-size system with 
periodic boundary conditions, the ground state is actu- 
ally a superposition of these two states, thus the value of 
the concurrence on the leg reduce to the half of the orig- 
inal value. In fact if we impose boundary condition in 
the same way as that in Ref. [3lJ , one of the two degen- 
erate states will be projected out, the staggered pattern 
of the leg concurrence appears and the value on the dim- 
mer leg is nearly 0.6 which is approaching to 0.7 in the 
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FIG. 3: (color online) The first derivative of the two-site en- 
tanglement on a rung as a function of 6. Lines of different 
system sizes cross at the same point (arctan(l/2), 0). 



rung singlet phase. In both the scalar chirality phase 
and dominant vector chirality phase, the concurrence of 
any pair of spin vanishes. However, at the cross-over re- 
gion between the dominant vector chirality and dominant 
collinear spin region, an unexpected concurrence on the 
diagonal pair appears and its maximal point (0 ~ 0.857r) 
is roughly the crossover point between the dominant vec- 
tor chirality region and dominant collinear spin region. 



IV. TWO-SITE ENTANGLEMENT OF THE 
RUNG AND THE SU(A) POINT 

In this section and the following section, we study the 
block-block entanglement of various blocks in this sys- 
tem. The block-block entanglement is the von Neumann 
entropy E v of a block of spin in the system. Precisely, it 
is calculated as: 



E V (A) = -tr (pa log 2 pa) 



(3) 



where A is a set of sites and pa is the corresponding 
reduced density matrix. If the whole system is in a pure 
state, then 



E V (A) = E V (B) = -tr [p B log 2 p B ) 



(4) 



where B is the rest part of the system. Then E v (A) or 
E V (B) describes how much the block A and the rest of 
the system are entangled. 

Compared to the concurrrence, the block-block entan- 
glement can apply to systems with much higher degrees 
of freedom, however, it is only meaningful when the con- 
cerning state is a pure state. In our calculation of finite- 
size ladders, it is found that the ground state is always 
non-degenerate in the region — 0.407T < 9 < 0.95-7T. Con- 
sidering the SU(2) symmetry of the Hamiltonian, the 
ground state's total spin is also zero in this region. 



The term two-site entanglement of a rung means that 
the von Neumann entropy is calculated from the reduced 
density matrix of two spins on the same rung. In Fig. 
[3] we show the first derivative of the entanglement as a 
function of 9 for different system size. From the figure, we 
observe that the first derivative of the entanglement on 
a rung reaches zero exactly at 9 = arctan (1/2) ~ 0.148-7T 
which is the QPT point between the staggered dim- 
mer phase and the scalar chirality phase. We find that 
this result is independent of system size. Recently, it 
was pointed out that, at this QPT point the system 
restores 5C/(4) symmetry p9l|. Precisely speaking, at 
9 C = arctan(l/2), the Hamiltonian commutes with the 
following operator [2jJ: 



T — J]] "5*1,1 • $2, j. 



(5) 



We will show that the expectation value of T is maximal 
or minimal exactly at 9 = 9 C due to the above symmetry. 

As discussed above, we can assume the ground state 
l^o) is non-degenerate around 9 Cl which implies \ipo) is 
also the eigenstate the of T at 9 C , thus we have T \ipo) = 
Xt \tj>o) i n which X t is some real number. Then the first 
derivative of (T) with 9 at 8 C is 



d0 



<Vo|T|Vo> = 



{ —\T\^) + (^\T\ — ) 



(6) 



= A t (— I*,) + A. <*,|_> + 



= 0. 



(7) 



Therefore, the expectation value of T reaches local 
maximum or minimum at 9 C . Since the ground state 
has j ky ) = (0, 0), the system is invariant under trans- 
lation along the leg, so (Sij ■ Sa.i) = wi,j ■ §2,j) (for 
any two site i and j in the ladder) is either maximal or 
minimal at 9 C . 

Next we show that there is one-to-one correspondence 
between the (Si^ ■ S2, i) and the two-site entanglement 
on the rung in the vicinity of 9 C . Let pij be the reduced 
density matrix of spins of sites i and j of the ground state. 
In the basis {|||) , |J.T) , |TI) , ITT)}, Pij has the following 
form due to the U(l) symmetry of the ground state in 
the concerning region(0.l7r < 9 < 0.27r): 



Pij 



/ u+ 

wt z* 

z w 2 

V u~ 



(8) 



Moreover, in the vicinity of 9 = 6 C , the ground state is 
unique and its total spin S = 0, which implies the ground 
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(b) 




FIG. 4: Four choices of increasing the block size: (a) is the 
single block, (b) is the stripe block, (c) is the zigzag block, 
(d) is one-leg block. 



state is also invariant under any rotations. Particularly, 
Pij is invariant under the rotation around the x axis: 

[at + aJ, Pij ] =0. (9) 

From Eq. @, and the condition trp = 1, we have 



l + 2z 



l-2z 



Wi = w 2 



(10) 



3z 

T' 



Si 4 -S 2A ) = ~, (11) 

Then the block-block entanglment on the rung E r is 

E r = — 3u + log 2 u + — (wi — z) log 2 (uii — z). (12) 

From the above equations it is clear that this extremal 
behavior of two-site entanglement on the rung is directly 
related to the SU (4) symmetry. 



V. SCALING BEHAVIOR OF THE 
BLOCK-BLOCK ENTANGLEMENT 

In this section, we focus on the scaling behavior of the 
block-block entanglement, namely how the block-block 
entanglement behaves as the block size changes. Unlike 
the case in the one dimensional chain, the ladder geom- 
etry has provided us many choices of how to select the 
block's shape and how to increase the block size. As 
shown in Fig. 0] we choose four different ways to increase 
the block size. 

First, we notice that in the ferromagnetic phase, the 
value of the block-block entanglement depends only on 
the size of the block. It is independent of the block ge- 
ometry. This is because any two sites in this state is 
essentially equivalent as we have already seen in Sec. IIIII 
In fact an explicit expression of block-block entanglement 
as a function of block size I can be obtained in the fol- 
lowing calculation. 

The ferromagnetic S^ ot = J2 i S? = state ipFM(S^ ot = 
0) could be obtained by applying lowering operator on 
the ip FM (Sf ot = N/2) N/2 times. 



'JFM 



(0)} 



(Ef =1 



MY- 



ITTT •■■?>, 



(13) 



JV 



where S t is the spin lowering operator at site 
i and N is the total number of sites. In the 
basis {|||) , ||t) > IU) ! ITT)}) a ll the coefficien ts of 
i>FM (Stot = 0) are the same: y/{N/2)\{N/2)\/N\. The 
matrix element pi{p,q) could be obtained explicitly: 



Pl(p,q) = I (^-p.)'(^+P')'W 




if Pz = q z 

if p z ^ q z 



(14) 



In the above expression, pi is reduced density matrix 
of the block consisting of I spins(we assume I < N/2), p 
and q are the column and row index of p;, p z and q z are 
the corresponding S% ot number. 

After diagonalizing this matrix, there are only I + 1 
nonzero eigenvalues X Pz withp z = —1/2,— Z/2+1, ... ,1/2. 
Then the block-block entanglement can be obtained as: 



(0!(^)! (#)!(#)! 



Pz 



E V {1) 



1/2 

E 



- V lo S2 



(15) 



When I and N are large, the summation in Eq. i|15|) 
can be replaced by an integral, and the function X Pz can 
be approximated by the Gaussian Distribution. Thus we 
can approximate E v l by the following expression: 



E v (l) 



log 2 



1 



N -I 



\ i0 ^ (y)' 



(16) 



which suggests that E v (l) diverges logarithmically as the 
size I increases. 

Secondly, in Fig. [ST a), in the most region of rung 
singlet phase, the block-block entanglement converges 
to some finite value quickly, while inEfb)-(d), it is al- 
most proportional to the block size. Since the ground 
state is approximately the rung singlet product state, 
the block-block entanglement is proportional to the num- 
ber of bonds cross the boundary of the block such as in 
case (c) and (d). In case (b), the situation is different 
since the number of bonds cross the boundary of block 
is finite, thus the short-ranged correlations between the 
rungs plays important role in (b) and explains this pro- 
portionality. 

Next, in Fig. EJb)-(d), we find that some local ex- 
treme points of the block-block entanglement may be the 
QPT points. Previous studies[2^,|2^| suggested that the 
transition point between the rung singlet and staggered 





FIG. 5: (color online) The block-block entanglement of different block size I as a function of 9 in N = 12 x 2 ladder, 
geometry of the blocks in (a)(b)(c)(d) are specified in Fig. |%[a)(b)(c)(d). 



The 



dimmer to be 0.067T ~ 0.087T which is quite near the lo- 
cal maximum point 0.07-7T in (b), local minimum 0.057T in 
(c) and (d). The transition point between the staggered 
dimmer and scalar chirality phase is exactly 0.1487T which 
is also near one minimal point 0.147T in (b). As for the 
crossover point between the dominant vector chirality re- 
gion and the dominant collinear spin region, (b)-(d) all 
suggest the value ~ 0.85-7T which is coincidence with the 
value obtained in previous works pjjj. 

In general, the scaling behavior of the above four 
choices of blocks could be categorized into two kinds, Fig. 
02a) and (b)-(d). In (a), the number of ladder bonds 
across the boundary between two blocks is a finite value 
which equals to 4 independent of the block size, while 
in (b)-(d), the number is proportional to the block size 
/. In the latter case, the short-ranged correlation across 
the boundary bonds is main contribution to the block- 
block entanglement, thus the block-block entanglement 
is always proportional to the size of the block as we have 



seen in Fig. 0] In the former case, the main contri- 
bution to the block-block entanglement comes from the 
long-range correlation between the sites in the block and 
the sites outside the blocks. It is expected that around 
the QPT point, in the former case the scaling behavior 
changes abruptly, e.g. from convergence to finite value to 
divergence logarithmically, while in the latter case, there 
may exist extremal point of the block-block entanglement 
which is an indication of QPT. 



VI. SUMMARY AND ACKNOWLEDGMENT 

In summary, we have studied the concurrence and 
the block-block entanglement in the ground state of the 
5=1/2 spin ladder with ring exchange. In both the 
rung singlet and staggered dimer phase, the behaviors 
of the ground-state concurrence are consistent with the 
corresponding dominant configurations. The extremal 
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point of the two-site block-block entanglement coincide 
with the QPT point due to the SU(A) symmetry, and 
such a symmetry is obviously independent of the system 
size. We have also investigated the scaling behavior of 
the block-block entanglement for different block geom- 
etry and block size. We have identified three kinds of 
typical scaling behavior in this model, namely increasing 
the size of block, the block-block entanglement (a) con- 
verges to some finite value, (b) diverges logarithmically 



with size, (c) diverges proportional with the size. How- 
ever, as we can see that, there's no signature of the QPT 
between the scalar chirality phase and dominant vector 
chirality phase. 
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for Research from the Research Grants Council of 
HKSAR, China (Project CUHK N_CUHK204/05 and 
HKU_3/05C). 
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